      program lesomabin
!        le duas matrizes em arquivos binarios soma e imprime arquivo binario
      implicit none
      integer, parameter:: nlinhas=900, ncolunas=900
      integer:: i,j
      real:: M1(nlinhas,ncolunas)
      real:: M2(nlinhas,ncolunas)


      open(1,file='bin1.bin',status="unknown",form="unformatted",access="direct",recl=4*nlinhas*ncolunas)
      read(1,rec=1)((M1(i,j),i=1,nlinhas),j=1,ncolunas)
      close(1) 

      open(1,file='bin2.bin',status="unknown",form="unformatted",access="direct",recl=4*nlinhas*ncolunas)
      read(1,rec=1)((M2(i,j),i=1,nlinhas),j=1,ncolunas)
      close(1)  


      do i=1,nlinhas
	  do j=1,ncolunas
	      M1(i,j)=M1(i,j)+M2(i,j)
	  enddo
      enddo

      open(1,file='bin3.bin',status="unknown",form="unformatted",access="direct",recl=4*nlinhas*ncolunas)
      write(1,rec=1)((M1(i,j),i=1,nlinhas),j=1,ncolunas)
      close(1)  


      stop
      end